

function domain = CalcDomainBlock(mesh)
nDim = mesh.nDim;
domain = zeros(1,nDim*2);

for iDim = 1:nDim
    domain(iDim) = min(mesh.coordinate(:,iDim));
    domain(iDim+nDim) = max(mesh.coordinate(:,iDim));
end